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Abstract 

We study a class of models with long-range repulsive interactions of the generalized Coulomb form V(r) ~ l/r Q . 
We show that decreasing the interaction exponent in the regime a < d dramatically depresses the charge order- 
ing temperature T c in any dimension d > 2, reflecting the strong geometric frustration produced by long-range 
interactions. A nearly frozen Coulomb liquid then survives in a broad pseudogap phase found at T > T c , which 
is characterized by an unusual temperature dependence of all quantities. In contrast, the leading critical behavior 
very close to the charge-ordering temperature remains identical as in models with short-range interactions. 
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1. Introduction 

The phenomenon of screening has long been known 
in presence of long-range Coulomb interactions, and it 
is generally expected to render the observed behavior 
very similar to that of systems with short-range inter- 
actions. Indeed, recent computational work has inves- 
tigated the critical behavior close to charge ordering in 
lattice Coulomb systems [1], suggesting the Ising uni- 
versality class. In addition, analytical arguments have 
been presented [2] supporting this view, for a broad 
class of lattice models with the generalized Coulomb in- 
teraction of the form V(r) ~ l/r a in d > 2 dimensions. 
We should mention, however, that long-range interac- 
tions are generally expected to produce mean-field like 
critical behavior for d < a < d* , while for d < d* < a 
one expect short-range critical behavior. The nontriv- 
ial effect of long-range interactions upon the critical 
behavior is possible only for a > d, because in this 
regime the interaction is "integrable" (no neutralizing 
background is needed) , and screening becomes inoper- 
ative. 

Should one expect any interesting or novel physics 
in Coulomb- like models (a < d), in comparison to the 
short-range situation? Our work confirms that these 



models indeed feature conventional critical behavior in 
the narrow critical region T ~ T c . We show, however, 
that a striking new behavior is uncovered in a broader 
temperature interval, reflecting strong geometric frus- 
tration inherent to such long-range interactions. First, 
we find a dramatic decrease of melting temperature of 
lattice Coulomb g result of the level of frustra- 

tion in the system. This can be easily understood by 
noting that our lattice Coulomb gas maps into an an- 
tiferromagnetic Ising model with long-range interac- 
tions. Here, all the spin tend to anti-align with all other 
spins, but this cannot be achieved for very long inter- 
actions (a < d), resulting in very low melting temper- 
ature. Indeed, for half-filled Coulomb systems (a = 1), 
the melting temperature T c is one order of magnitude 
smaller [3] than the generalized Coulomb energy E c — 
e 2 /a , where a is the lattice constant. Continuum mod- 
els [4-6] show an even more dramatic behavior, with 
the melting temperature being as much as two orders 
of magnitude smaller then the Coulomb energy. This 
striking behavior, although well documented in several 
model studies, is not widely appreciated or understood 
in simple terms. 

The second robust feature of these models, which has 
only recently been discovered [7], is the emergence of 
the "pseudogap phase". This pseudogap phase is a spe- 
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cific feature of long-range type interaction with a < d, 
and is observed in a broad temperature range T c < 
T < T* (see Fig. 1), where T* is the pseudogap tem- 
perature where the gap in the single particle density of 
state (DOS) starts to open. We show that the physical 
Coulomb interaction (a = 1) lays deep in the regime 
of very long range interaction a — > 0, the regime where 
our analytical "Extended Dynamical Mean-Field The- 
ory" (EDMFT) approach becomes is asymptotically 
exact [7]. This observation explains the surprising ac- 
curacy of this analytical scheme both when applied 
to clean models [7], and in previous applications to 
Coulomb glasses [8]. 



2. Model and EDMFT Approach 

Our model is described by the Hamiltonian of par- 
ticles living on a half-filled hypercube lattice (lattice 
spacing a) that interact via Coulomb-like interaction 

H = \Y. v (Rij)(m - (n»K - ("))■ (i) 

ij 

Here Rij is the distance between lattice sites i and j 
expressed in the units of the lattice spacing. In general, 
our unit of energy the nearest-neighbor repulsion E c = 
V(l) . We focus on the half-filled system < n >= 1/2 
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Fig. 1. (Color online) Phase diagram of the half-filled clas- 
sical D — 3 lattice model with interactions V(R) — R~~ a . 
The charge ordering temperature T c (a) ~ a, as obtained from 
EDMFT theory (full line) and Monte-Carlo simulations (open 
symbols) . The pseudogap temperature T* (dashed line) remains 
finite as a —> 0; a broad pseudogap phase emerges at a < d. 
We also show T^ R ?a 1 for the same model with short-range 
interactions (dotted line). The inset shows corresponding sin- 
gle particle density of state (DOS) for different temperatures. 
The EDMFT results (dashed line) show excellent agreement 
with Monte-Carlo simulation (solid line) above the melting 
temperature T c . [7,8] 
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Fig. 2. (Color online) Scaling behavior of the spectral den- 
sity v(e) of plasmon modes on a 3D hypercubic lattice. The 
low-energy branch of the spectrum, which describes the "sheer" 
charge fluctuations close to the ordering wavevector, features 
an a-dependont energy scale, one that sets the melting temper- 
ature T c ~ a. The scaling becomes exact for very long inter- 
action range a <^ 1, but is approximately valid for all a < 1, 
explaining the simple a-dcpcndcncc of all quantities. 

that maps the problem into "antiferromagnetic-like" 
Ising model. In addition to numerically-exact (classi- 
cal) Monte-Carlo simulations, we also utilize an ana- 
lytical approach, called so-called Extended Dynamical 
Mean Field Theory (EDMFT) [9,10] which is expected 
[7] to be accurate for very long interaction (a <g 1), 
where the effective coordination number becomes very 
large. The inter-site charge correlations are included 
in this approach, and are treated on the same foot- 
ing as the local ones [10]. Here, the lattice problem is 
mapped into an effective impurity problem embedded 
in a self-consistently determined bath containing both 
fermionic and bosonic excitations [9,10]. In our case, 
this bosonic bath describes the plasmon excitations in- 
duced by the inter-site charge correlations. 

In this paper, we focus on the classical limit, where 
the origin of the pseudogap phase is most pronounced 
and the accuracy of the EDMFT can be tested directly 
with simple classical Monte Carlo simulations. Here, 
the density-density correlation function becomes 

X (k) = (4 + D + f5V k r\ (2) 
and the self-consistent condition 

i= J deiy(e)(4 + D + /3e)-\ (3) 

where we define the classical plasmon-mode spectral 
density (Fig. 2) u{e) = ~}2 k 8(e — 14). By solving nu- 
merically the self-consistent condition in equation (3) 
using the spectral density v{e), the T dependence of 
dimensionless parameter D(T) is found to be power- 
law dependent (Fig. 3). The analytic solution can be 
derived at the limit of a d, spectral density becomes 
v(e) ^ S(e — so), where e = Vq {Q = {it, it) in 2D, 
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Fig. 3. (Color online) The power law behavior of D(T) for 
different interaction range a. In the limit of a — > 0, we find 
parameter D(T) ~ T -1 , both in d=2 and d=3. 

and Q = (tt, tt, tt) in 3D) is the minimum energy in mo- 
mentum space (center of Brillouin zone). By solving 
equation (2) at this limit, we find the analytic solution 
D(T) = T~ % Vq. The result qualitatively holds all the 
way to a ~ d as shown in Fig. 3. 

The charge ordering temperature T c (a) is defined by 
the x(^) — > oo at the corresponding ordering wave vec- 
tor k = Q; in d = 3 the system forms a BCC structure, 
and a checkerboard patter in d — 2. The depression of 
melting temperature as a function of interaction pa- 
rameter a < 1 can be understood in simple way by 
noting that the spectral density v(e) has a simple scal- 
ing form 

v(e) = a~ 1 ?{(e - e )/a}, (4) 

as shown in Fig. 2. The sharp peak of (e — so) below 
characteristic energy scale e*(a) ~ a, physically cor- 
respond to "sheer" plasmon modes with wave vector 
k ~ Q. [5,11] This energy scale plays a role of an effec- 
tive Debye temperature. This tells us why the ordering 
temperature decreases 



T c (a) = a J deu(e)/e ~ e*(a) 



0.1a, 



(5) 



in agreement with the estimation of Lindermann cri- 
terion applied to the sheer mode. At the lowest energy 
near so the dispersion relation assumes the standard 
form v(e) ~ £ ( d_2 )/ 2 i same as for short-range interac- 
tions. This can also be understood from the form of 
the potential surface in momentum space, as shown in 
Fig. 4. The potential surface is flatter for smaller a 
and the values are close to the minimum energy level 
at Q = (tt, tt). This picture clearly explains the origin 
of the sharp peak in the spectral density u(e). 

In the classical limit, the single particle density of 
state (DOS) is nothing more than the probability dis- 
tribution function of electrostatic potentials computed 
on different sites in the lattice 



p(u,,T) = (^2S(u -52v ijnj )J . (6) 
Within the EDMFT approach it assumes the form 
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The Coulomb pseudogap ("plasma dip") E gap = 
D/P , which is the distance between the Gaussian 
peaks that start to open at D — E gap /4: ■ As shown 
in Fig.l, both E gap and T* remain finite for a <C 1, 
since D(T) ?s j3 in this limit (Fig. 3). This results, 
for long range interactions, in the formation of a wide 
pseudogap phase, since the melting temperature de- 
creases as T c (a) ~ 0.1a while pseudogap temperature 
remains finite (T* — > 1/4) as a -> 0. We define T* 
based on the change of curvature in the DOS at u = 0, 
i.e. the vanishing of the second derivative of the DOS 
d 2 p{uj)/du = -> D = 8, and finally use D(T*) = 8 
from the power- law dependent of D(T) (Fig. 3). 



3. Long Range Interaction: Universality Class 

The fact that models with long-range interactions 
with a < d belong to the same universality class as the 
short-range Ising model can be clearly seen from care- 
fully examining the form of the interaction in momen- 
tum space. To determine its precise form in quantita- 
tive detail one must use Ewald summation [12-14] and 
split the appropriate lattice into a short-range part 
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and long-range part 
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where F(a, x) is the incomplete Gamma function, and 
G is reciprocal lattice unit vector. 

The resulting dispersion is shown, for different a 
in Fig. 4. At the center of the Brillouin zone, the 
long-range part for different interaction a is equally 
parabolic for all k vectors and only start to drop 
sharply when approaching the vicinity of the corners of 
momentum space Z*. The singularity points near the 
edges of the Brillouin zone where fej = 2ttZ, diverges 
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to Vk — » oo, leading to positively diverging energy for 
corresponding (longitudinal) fluctuation modes. This 
mode has very small statistical weight. It is related to 
the tail on the spectral density (Fig. 2), and it does 
not contribute in the critical behavior of the system. 
The critical behavior, in contrast, is determine by the 
dispersion of "sheer" modes in the vicinity of the order- 
ing wavevector, a quantity that remains of the same 
qualitative form as for short-range interactions [2]. 




Fig. 4. (Color online) (a) The dispersion of the long-range 
part of Coulomb potential in d — 2 (Jo — —1 [re- 
verse potential surface]) according to equation (8,9) for 
different interaction range. The lowest energy state near 
Q — (7r, 7r) feature parabolic like surface for interactions a 
equal to 1 (red/flattest), 2(green/middlc),3(blue/curviest) and 
only change rapidly near the edges of the first Brillouin zone. 
The edges arc related to the tail on the plasmon spectral den- 
sity of Fig. 2 and the parabolic surface related to the peak on 
the plasmon spectral density. At a > d, the long-range dis- 
persion expansion docs not have rapid change near the edges 
and is more like short-range dispersion, (b) Short-range dis- 
persion expansion does not have singularity near the edges for 
any interaction range. 



4. Conclusions 

The role of long-range interaction in suppressing crit- 
ical temperature T c is an important feature of systems 
with long-range interactions, behavior which directly 
reflects the generic presence of strong geometric frus- 
tration. While these effects do not modify the univer- 
sality class of the narrow critical region, novel behavior 
is found well above the melting temperature. Here we 
uncovered the emergence of a broad pseudogap phase, 
a feature unique to long-range interactions; this regime 
is very accurately described by our analytical EDMFT 
theory, and we find excellent agreement with (numeri- 
cally exact) Monte Carlo simulations. 
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